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ABSTRACT 

We present Infrared Array Camera (IRAC) data and source catalogs from the Spitzer Space Tele- 
scope Extragalactic First Look Survey. The data were taken in four broad bands centered at nominal 
wavelengths of 3.6, 4.5, 5.8 and 8.0/zm. A set of mosaics and catalogs have been produced which are 
« 80% complete and « 99% reliable to their chosen flux density limits. The main field survey covers 
3.8deg^, and has flux density limits of 20/xJy, 25/iJy, 100/xJy and 100/iJy at wavelengths of 3.6, 4.5, 
5.8 and 8.0/im, respectively. The deeper "verification" survey covers 0.25deg^ with limits of 10/iJy, 
lO^Jy, 30^Jy and 30/iJy, respectively. We also include deep data in the ELAIS-Nl field which covers 
0.041deg^ with limits of 4/iJy, 3/iJy, lO/iJy and 10/iJy, respectively, but with only two wavelength 
coverage at a given sky position. The final bandmerged catalogs contain 103193 objects in the main 
field, 12224 in the verification field and 5239 in ELAIS-Nl. Flux densities of high signal-to-noise 
objects are accurate to about 10%, and the residual systematic error in the absolute flux density scale 
is ~ 2 — 3%. We have successfully extracted sources at source densities as high as 100000 deg^^ in 
our deepest 3.6 and 4.5/j,m data. The mosaics and source catalogs will be made available through the 
Spitzer Science Center archive and the Infrared Science Archive. 
Subject headings: surveys - catalogs - infrared radiation 



1. INTRODUCTION 

The extragalactic portion of the Spitzer Space Tele- 
scope (Werner et al. 2004) First Look Survey (here- 
after XFLS) was one of the first observations made with 
Spitzer after the completion of Science Verification at the 
end of November 2003. The aim of this 67hr survey was 
to characterize the extragalactic source populations ob- 
served with Spitzer down to sub-milliJansky levels in the 
mid- infrared Observations covering the survey areas 
were made with both the Infrared Array Camera (IRAC; 
Fazio et al. 2004a) and the Multiband Imaging Photome- 
ter for Spitzer (MIPS). This paper discusses the IRAC 
observations, the MIPS observations will be discussed in 
future papers (Fadda et al. in preparation; Frayer et al. 
in preparation). IRAC takes images in four broad bands 
(termed channels). Channel 1 is centered at a nominal 
wavelength of 3.6/im, channel 2 at 4.5/im, channel 3 at 
5.8/im and channel 4 at 8.0/im Light enters the in- 
strument through two apertures, one for channels 1 and 
3, and one for channels 2 and 4. The centers of the two 
apertures are separated by about 6.'5 on the sky. Beam- 
splitters split the short and long wavelength light in each 
aperture. The light in channels 1 and 2 is detected by 
two 256 X 256 InSb arrays, and that in channels 3 and 
4 is detected by two 256 x 256 SiAs arrays. All arrays 
have the same pixel size, corresponding to 1'.'22 per pixel, 

^ Spitzer Science Center, Caltech, Mail Code 220-6, Pasadena, 
CA 91125; mlacy@ipac.caltech.edu 

^ Astronomy Program, School of Earth and Environmental Sci- 
ences, Seoul National University, Shillim-dong, Kwanak-gu, Seoul, 
S. Korea 2-880-9010 

^ See the First L ook Survey website at 
|http: / /ssc. spitzer. caltech.edu/fls/ 

The actual central wavelengths are slightly different and de- 
pend on the source spectrum, see Fazio et al. (2004) for details. 



giving a 5'.2 x 5'.2 field of view. The pipehne-processed 
IRAC data from the XFLS were publicly released at the 
end of April 2004. In this paper, we describe the data 
analysis which we performed to improve the data quality 
beyond the scope of the standard pipeline processing, in- 
cluding a discussion of the removal of artifacts from the 
data. We present IRAC source catalogs for the XFLS, 
and discuss good observing practices for future surveys 
of this nature. 

2. DATA COLLECTION 

Nine Astronomical Observation Requests (AORs) 
cover the main field survey area of 3.8deg^ centered on 
R.A.(J2000) 17'^18'"00^ Dec.(J2000) -t-59°30'00" in a 
3x3 grid in array (row, column) coordinates. Each AOR 
was an approximately 8x8 map with 277 offsets. The 
small 5-point Gaussian dither pattern with 12s frametime 
was used. This pattern has a mean offset of 28 . Data 
in the verification area (a 0.25 deg^ area within the main 
field centered on R.A.(J2000) 17'^17"^00^ Dec.(J2000) 
-|-59°45 00 ) was taken using three AORs with 12s fram- 
etimes with the same dither and mapping strategy as the 
main field and three much deeper AORs with 30s frame- 
times using the first sixteen points from the small cycling 
dither pattern, which has a mean offset of 13 . These 
AORs are all contained in Spitzer program identification 
number (PID) 26 (principal investigator (P.I.) T. Soifer). 

In addition, a second field in the ELAIS-Nl region, cen- 
tered on R.A.(J2000) 16''09'°20% Dec.(2000) -h54°57'00" 
was observed as part of a study of source confusion as a 
test field for the Great Observatories Origins Deep Sur- 
vey (GOODS), and this has also been included in the 
XFLS (PID 196, P.I. M. Dickinson). These data con- 
sisted of a two position map separated by 312 , with the 
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first 36 dithers from the medium cycling pattern and a 
200s frametime. The map was oriented to avoid overlap 
between the channel 1/3 and channel 2/4 fields of view, 
thus, unlike the main and verification fields, a typical sky 
position in the image only has two channel coverage. 

Details of the XFLS AORs, along with their AOR iden- 
tification numbers (AORIDs) are given in Table 1. The 
AORs can be downloaded using the Spot tool (avail- 
able from http://ssc.spitzcr.caltech.edu/propkit/spot/), 
selecting the "View Program" option and specifying the 
appropriate PID. 

3. DATA ANALYSIS 

3.1. Pipeline Processing 

The main and verification field data were run through 
the S10.5 version of the Spitzer Science Center (SSC) 
pipeline, described in the IRAC Data Handbook^. Those 
for the ELAIS-Nl field used the S9.5 version of the 
pipeline. Between S9.5 and S10.5 the treatment of the 
darks was improved, but the long frames in the ELAIS- 
Nl data are background dominated, so the improvement 
between S9.5 and S10.5 was neghgible for these data. 
For each IRAC frame (termed a Data Collection Event 
(DCE)) the IRAC pipeline produces processed images, 
the Basic Calibrated Data (BCD), and corresponding 
masks (the DCE masks, or Dmasks) which we used as 
the starting point for our analysis. The masks flag po- 
tential problems with the data such as saturated pixels, 
strong radiation hits or corrupted/missing data. 

3.2. Post-pipeline Processing 

The IRAC data contain a number of artifacts which 
are not yet removed well by the SSC pipeline. Some of 
these were removed in post-pipeline processing steps. 

Bright stars and cosmic rays systematically decrease 
the bias level in the columns in which they lie in chan- 
nels 1 and 2 (Figure 1). This "column pulldown" was 
removed using an algorithm developed and implemented 
by L. Moustakas and D. Stern (personal communication). 
This masks objects in the frame to produce a background 
image, then searches for discrepant columns, correcting 
them by adding a constant to bring the mean background 
value in the column up to that of its neighbors. Multi- 
plexer bleed ( "muxbleed" ) in channels 1 and 2 from stars 
and cosmic rays with peak fluxes in the linear regime 
was partly corrected by the SSC pipeline, but that from 
the brightest, saturated stars was not. To correct these, 
we searched for bright stars in the frame, then split the 
data into the four amplifier outputs (every fourth column 
of data passes through the same output amplifier). For 
each of the four outputs, we fit the rows within five pixels 
of the row containing the star with a straight line using 
columns outside ±30 pixels of the star, subtracted the fit, 
and reassembled the frame. In channels 3 and 4, band- 
ing effects (predominately optical in origin, produced by 
scattering of light within the detector) were present (Fig- 
ure 2). These resulted in the rows and columns contain- 
ing the bright stars to be elevated above the background. 
The banding artifacts were removed in a similar fash- 
ion to the pulldown and muxbleed in channels 1 and 2, 
though in this case the outputs were not split. Muxbleed 

Available from the SSC website 

(ssc.spitzer.caltech.edu/irac/dh/) 



from extremely bright stars in channels 1 and 2 can re- 
sult in a large fraction of the array being offset from the 
mean level. This was corrected by searching for saturated 
stars, then averaging the rows above and below the star 
in each of the four outputs, fitting a line to the mean of 
each row and subtracting the fit. Further details on and 
examples of these artifacts may be found in the IRAC 
Data Handbook. 

The IRAC data have a number of problems which af- 
fect the background level across frames. The "first frame 
effect" - the variation in the dark which depends on the 
length of time since the previous exposure - is partly cor- 
rected in the pipeline, but some residual variation in the 
dark remains, particularly in channel 3. There are two 
reasons for this. First, the first frame calibrations (taken 
prior to launch), did not cover a sufficiently large time 
range to calibrate the longest-term variations, and sec- 
ond, they were taken prior to the decision to anneal the 
IRAC arrays regularly to remove long-term latents. To 
remove the dark variations to first order, a plane was fit- 
ted to produce a flat image with a constant background. 
This dealt with much of the uncorrected first frame ef- 
fect. However, background variations remained from 
three sources which produce background variations on 
scales smaller than the array. These were: residual first- 
frame/array relaxation effects, long-term latents (partic- 
ularly noticeable in channel 1), and some real positional 
variations in the background level due to high Galactic 
latitude dust emission in channel 4. We therefore sub- 
tracted a "delta dark" from each frame. In the data taken 
with the 12s frametimes, channels 1, 2 and 4 were ade- 
quately corrected (in the sense that the variations in the 
background level were reduced to be similar in magnitude 
to the noise level) by subtracting a ±1.5(T-clipped median 
of the whole pixel stack from all the frames. Channel 3 
required a running median of 15 of the ±1.5cr-clipped 
pixel stack to produce acceptable results. The 30s fram- 
etime data were more problematic. Channel 1 proved 
to be particularly challenging. The small dither pattern 
proved too small to effectively remove latent images from 
bright stars. In addition, the longer frametimes lead to 
a gradual build-up of latents over the whole array. A 
running median of 24 of the ibl.Scr-clipped pixel stack 
produced good results for all but the last two pointings 
of the AORs (which happened to contain bright stars). 
These frames had a 48-frame running median subtracted 
instead. In channel 3 the first 30s frame of the AOR 
was essentially uncorrectable and was excluded from the 
mosaics. For the rest, a 36-frame running median of the 
zbl.Scr-clipped pixel stack was found to be adequate. A 
straight median of the clipped frames, as used for the 12s 
data, was found to be adequate in channels 2 and 4. 

To mask scattered light from sources off the detector 
array^, a straylight masking program developed by R. 
Arendt and ML was used to place masks into the Dmask 
file. This program uses a list of bright stars from the 
Two Micron AU Sky Survey (2MASS; Cutri et al. 2003) 
to predict the positions of light scattered from the focal 
plane assembly (FPA) cover (the principal source of scat- 
tered light in channels 1 and 2), and scattering off the 
edge of the detector (the principal source of scattering 

^ For details of scattered light in IRAC see the 5pitzer O bservers' 
M anual |http://ssc.spitzer.caltech.edu/documents/som/| 
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AORID 



Date Obs 
(UT) 



Field Frametime 
(seconds) 



TABLE 1 

Spitzer OBSERVATIONS OF THE XFLS FIELDS 

Map offset* Mapping strategy^ Dither Pattern"^ 
(arcsec) 



3861504 

3861760 

3862016 

3862272 

3862528 

3862784 

3863040 

38633296 

3863552 

3866880 

7676928 

7677184 

3867136 

3867392 

3867648 

6006016 



2003- 
2003- 
2003- 
2003- 
2003- 
2003- 
2003- 
2003- 
2003- 
2003- 
2003- 
2003- 
2003- 
2003- 
2003- 
2003- 



12-01 
12-01 
12-02 
12-02 
12-03 
12-03 
12-03 
12-04 
12-04 
12-05 
12-05 
12-05 
12-05 
12-05 
12-06 
12-28 



mam 

main 

main 

main 

main 

main 

main 

main 

main 
verification 
verification 
verification 
verification 
verification 
verification 
ELAIS-Nl 



12 
12 
12 
12 
12 
12 
12 
12 
12 
12 
12 
12 
30 
30 
30 
200 



-2273.5,-2550.5 
-2412.0,-138.5 
-2412.0,2135.0 
-138.5,-2412.0 
-138.5,-138.5 
-138.5,2273.5 
2273.5,-2273.5 
2273.5,-138.5 
2135.0,2412.0 

0.0,0.0 
-830.0,-150.0 
830.0,150.0 

0,0 
830.0,150.0 
-830.0,-150.0 
201.65,0.0 



7x9, 
8x9, 
8x8, 
9x8, 
9x9, 
9x9, 
9x7, 
9x9, 
8 X 10 
4x7, 
2x6, 
2x6, 
4x5, 
2x5, 
2x5, 
1x2, 



(277 , 
(277", 
(277", 
(277", 
(277", 
(277", 
(277", 
(277", 
(277" 
(277", 
(277", 
(277", 
(277", 
(277", 
(277", 
(292", 



277 ) 
277") 
277") 
277") 
277") 
277") 
277") 
277") 
, 277") 
277") 
277") 
277") 
277") 
277") 
277") 
312") 



5-point Gaussian 
5-point Gaussian 
5-point Gaussian 
5-point Gaussian 
5-point Gaussian 
5-point Gaussian 
5-point Gaussian 
5-point Gaussian 
5-point Gaussian 
5-point Gaussian 
5-point Gaussian 
5-point Gaussian 
Points 1-16 from 
Points 1-16 from 
Points 1-16 from 
Points 1-36 from 



, small 
, small 
, small 
, small 
, small 
, small 
, small 
, small 
, small 
, small 
, small 
, small 

small cycling 
small cycling 
small cycling 
medium cycling 



* relative to the appropriate field center, in array coordinates. 

^ size of map (row X column), mapping steps in array coordinates. 

* these are standard IRAC dither patterns, described in the Spitzer Observers' Manuarhttp://ssc. spitzer. caltech.edu/documents/somj 




Sum 




Fig. 1. — A bright star in the channel 1 mosaic shown before (left) and after (right) the post- pipeline processing was applied. The 
star shows all the artifacts produced by a bright star in channels 1 and 2: column pulldown, the reduction in the data values in columns 
directly above and below the star, muxbleed (raising the level of pixels in rows either side of the star), and the offset in the background 
level produced by extremely strong muxbleed affecting all columns below the star, and also producing striations in the background level. 
The post-pipeline processing is able to remove most of the artifacts. 



in channels 3 and 4). These masks have been apphed to 
the final mosaics. 

Software for IRAC artifact correction and 
masking, including up to date versions of much 
of the code described above may be found 
on the Spi tzer contributed software website 
(http://ssc.spitzer.caltech.edu/archanaly/contributed/ 



browse.html). 

The SSC pointing refinement module (Masci, Makovoz 
& Moshir 2004) was run on all the channel 1 and chan- 
nel 2 BCDs. This module improves the pointing of the 
BCDs by both registering each BCD with respect to over- 
lapping BCDs ("relative refinement") and by registering 
the result with stars from the 2MASS catalog ("absolute 
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Fig. 2. — A bright star in the channel 3 mosaic shown before (left) and after (right) the post-pipeline processing was applied. The star 
shows the combination of optical banding and electronic artifacts produced by a bright star in channels 3 and 4, resulting in the raised 
values of the pixels in the rows and columns on which a bright star falls. In the channel 3 data the artifacts are more pronounced in the 
column than the row direction, but in channel 4 the stronger artifacts are in the row direction. See the IRAC Data Handbook for more 
details. 



refinement"). To apply these corrections to the channel 
3 and 4 BCDs, which have too few stars for a reliable ab- 
solute correction, we computed a weighted average of the 
channel 1 and channel 2 correction offsets, applied these 
new offsets to all channel 1,2,3 and 4 BCDs and finally 
computed the refined pointing for all the BCDs. The 
mosaics constructed using the refined pointings from the 
BCD are well-registered to the 2MASS frame (see Section 
5 below). 

3.3. Production of the Mosaics 

The SSC mosaicing software Mopex^ was used to mo- 
saic the individual BCD images. Mopex consists of a set 
of modules combined in a perl wrapper script. Which 
modules are used, and the parameters which control the 
modules, are determined by namelist files. The namelists 
used were similar to the one given in the IRAC Data 
Handbook. Standard linear interpolation was used for 
all the mosaics. Two modes of outlier rejection were em- 
ployed, the "dual outlier" rejection, which uses a com- 
bination of spatial and temporal criteria to identify out- 
liers, and a purely temporal ( "multiframe" ) technique 
which uses the pixel stack at a given image position to 
identify outliers. The multiframe technique is the most 
effective when the coverage is high (>4), but the dual out- 
lier method is better on the lower-coverage regions on the 
edge of the mosaics, or in regions masked by the Dmasks. 
The single frame outlier module, which relies purely on 
a spatial criterion to determine outliers, was not used. 
The diffuse cosmic rays in channels 3 and 4 required 
that the masks be grown around detected outlier pixels 
in the lower coverage data. This was achieved in prac- 
tice by setting the parameter RMJTHRESH to a low 
value (0.05) in the mosaicer channel 3 and 4 namelists for 
the main and verification fields. Coverage maps for each 
mosaic were output by the mosaicer, and were checked 
by eye for evidence of over-zealous outlier rejection. To 
help with dcblcnding sources the ELAIS-Nl data were 
mosaiced with a two-to-one pixel ratio to better sample 

Described in http://ssc.spitzer.caltech. edu / postbcd/docT] mo- 
saiker.pdf 



the data. The main field and verification strip retained 
their original sampling. The Mopex namelists used for 
the XFLS will be made available on the XFLS website 
so others can reproduce our results should they wish to 
do so. Some of the more important parameter choices 
are listed in Table 2. 

4. PRODUCTION OF THE SOURCE CATALOGS 

Single-band catalogs were produced from each mosaic 
image using Sextractor (Bertin & Arnouts 1996). The 
coverage maps produced by Mopex were used as weight 
images. We have optimized the photometry for the 
faintest objects. A background mesh size of 16 pixels 
with a filter width of three was used, and the "local" 
background (measured in a 24 pixel thick annulus) ap- 
plied. Large, extended objects are thus likely to have 
incorrect fiuxes in the catalog and should have their 
fiuxes remeasured from the mosaics. Four fixed aper- 
ture fluxes (aperture diameters 6'.'00, 9'.'26, 14'.'86 and 
24'.'40) plus an isophotal flux were measured. The IRAC 
point spread function has broad wings compared to typ- 
ical ground-based data, and aperture corrections are sig- 
nificant (see the IRAC Data Handbook for further de- 
tails). The largest aperture diameter (24'.'40) is the same 
diameter that the IRAC calibration stars are measured 
in. The fluxes of the calibration stars in this aperture are 
considered the total fluxes in the S9.5 and S10.5 pipefine 
processed data, so this aperture requires no correction. 
The smaller apertures were picked to match those used 
by the Sloan Digital Sky Survey (SDSS). For these, aper- 
ture corrections were applied using a lookup table de- 
rived from measurements of bright stars in the XFLS 
data, and checked for consistency with those derived by 
the IRAC instrument team and those in the IRAC Data 
Handbook (see Table 3). Although these corrections 
are not quite appropriate for the typical XFLS galaxy, 
which is slightly extended in IRAC images, our simula- 
tions showed that the application of this correction nev- 
ertheless significantly improved the flux densities. One 
of the aperture fluxes is taken as a "best" flux in the 
catalog. The "best" aperture was chosen by comparing 
the geometric mean radius of the isophote (r,„ ^JA/t: 
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TABLE 2 

Values of the more important mosaicer and Sextractor parameters 

Common to all: 



Program (module) 


parameter 


channel 1 


channel 2 


channel 3 


channel 4 


Mopcx (DETECT) 


Detect ion_M ax _Area 


9 


9 


9 


9 


Mopox (DETECT) 


Detect ion_Min_Arca 














Mopcx (DETECT) 


Detect ion_Threshold 


5 


5 


5 


5 


Mopcx (MOSAICINT) 


INTERP .METHOD 


1 


1 


1 


1 


Mopcx (MOSAICDUALOUTLIER) 


MIN.OUTL.IMAGE 


2 


2 


2 


2 


Mopcx (MOSAICDUALOUTLIER) 


MIN.OUTL.FRAC 


0.2 


0.2 


0.2 


0.2 


Mopcx (MOSAICOUTLIER) 


THRESH.OPTION 


2 


2 


2 


2 


Mopcx (MOSAICOUTLIER) 


BOTTOM.THRESHOLD 


4 


4 


4 


4 


Mopcx (MOSAICOUTLIER) 


TOP.THRESHOLD 


4 


4 


4 


4 


Mopcx (MOSAICOUTLIER) 


MIN.PIX.NUM 


3 


3 


3 


3 


Mopcx (MOSAICRMASK) 


MIN.COVERAGE 


3 


3 


3 


3 


Mopcx (MOSAICRMASK) 


MAX.COVERAGE 


100 


100 


100 


100 


Sextractor 


DETECT.MINAREA 


3 


3 


3 


3 


Sextractor 


DETECT.THRESH 


2 


2 


2 


2 


Sextractor 


ANALYSIS.THRESH 


2 


2 


2 


2 


Sextractor 


FILTER 


N 


N 


N 


N 


Sextractor 


DEBLEND.NTHRESH 


32 


32 


32 


32 


Sextractor 


DEBLEND.MINCONT 


0.0001 


0.0001 


0.0001 


0.0001 


Sextractor 


BACK.SIZE 


16 


16 


16 


16 


Sextractor 


BACK.FILTERSIZE 


3 


3 


3 


3 


Sextractor 


BACKPHOTO.TYPE 


LOCAL 


LOCAL 


LOCAL 


LOCAL 


Main field: 












Program 


parameter 


channel 1 


channel 2 


channel 3 


channel 4 


Mopex 


MOSAIC JIXEL.RATIOJC 


1 


1 


1 


1 


Mopex 


MOSAIC JIXEL.RATIO.Y 


1 


1 


1 


1 


Mopex (MOSAICRMASK) 


RM.THRESH 


0.8 


0.8 


0.05 


0.05 


ScxI.rartor 


SATTTR .LEVEL 


.'525 


in. '5 


2271 


878 


ScxUacLor 


CiAlN 


:iU9 


2(S0 


(i7 


IIJ.") 


Verification field: 












Program 


parameter 


channel 1 


channel 2 


channel 3 


channel 4 


Mopex 


MOSAIC JIXEL.RATIOJC 


1 


1 


1 


1 


Mopex 


MOSAIC JIXEL.RATIO.Y 


1 


1 


1 


1 


Mopex (MOSAICRMASK) 


RM.THRESH 


0.8 


0.8 


0.05 


0.05 


Sextractor 


SATUR.LEVEL 


126 


180 


883 


341 


Sextractor 


GAIN 


786 


721 


67 


172 


ELAIS-Nl field: 












Program 


parameter 


channel 1 


channel 2 


channel 3 


channel 4 


Mopex 


MOSAIC.PIXEL.RATIOJC 


2 


2 


2 


2 


Mopcx 


MOSAIC.PIXEL.RATIO.Y 


2 


2 


2 


2 


Mopcx (MOSAICRMASK) 


RM.THRESH 


0.8 


0.8 


0.8 


0.8 


Sextractor 


SATUR.LEVEL 


17.4 


24.9 


122 


195 


Sextractor 


GAIN 


5679 


5210 


1244 


876 



thus reduces noise and confusion in the flux density mea- 
surement. 

Flag images based on the results of the muxbleed, pull- 
down and banding correctors, latent image masking in 

the final mosaic (applied by hand) and haloes of bright 
stars (using sizes predicted from 2MASS K-magnitudes) 
were applied to the Sextractor catalogs. Objects whose 
isophotes cross non-zero pixels in the flag image have 
a value in the image flag field derived from OR-ing to- 
gether the fiag image pixels which fall within the object 
isophotes. Flag values and their meaning, for both the 
flag values output by Sextractor and those added by ap- 
plying the flag image, are listed in Table 4. 

The four IRAC single-band catalogs were also merged 
into a single. 4-band catalogue. The matching procedure 
began by going through each source in the band 1 catalog 
and searching for a match within a radius of 1.5 pixels 
in bands 2-4. Sources in band 2 unmatched to band 1 



TABLE 3 
Aperture corrections 

Aperture Correction applied 

channel 1 channel 2 channel 3 channel 4 
Woo 1.167 1.213 1.237 1.466 
9'.'26 1.091 1.117 1.100 1.165 
14'.'86 1.042 1.048 1.042 1.066 
24^'4 1.000 1.000 1.000 1.000 



where A is the isophotal area) to each fixed aperture di- 
ameter. This radius is compared to the radii of the four 
apertures 'ri, r2,r3 and r^. If r„i < l.lri then the "best" 
aperture radius, r^, is set to ri, if l.lri < < l.lr2 
then Tb = r2 etc. If > l.lr4 then the isophotal flux is 
used as the "best" flux. Although somewhat arbitrary, 
this procedure ensures that an aperture appropriate to 
the isophotal size of the source is selected as "best" , and 
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TABLE 4 
Image flag values 



Bit meaning 



set by 




1 
2 
3 
4 
5 
6 
7 
8 
9 
10 



More than 10% area affected by bad pixels 

Object originally blended with another one 

At least one object pixel is saturated 

Object is truncated at image boundary 

Object's aperture data are incomplete or corrupted 



Sextractor 
Scxtractor 
Sextractor 
Sextractor 
Sextractor 



but the field is sufficiently deep that IRAC and ACS 
data from the Great Observatories Origins Deep Sur- 
vey (GOODS) in the Hubble Deep Field North (HDFN) 
-eould be used instead. Ten galaxies close to the flux 



Object's isophotal data are incomplete or corrupted Sextractor 



A memory overflow occurred during deblending 
A memory overflow occurred during extraction 
Region aflfected by muxbleed, pulldown or banding 
Object in halo of bright star 
Object contaminated by latent image 



Sextra<;tor 
Sextractor 



sources were then matched to bands 3 and 4, and finally 
sources in band 3 unmatched to bands 1 and 2 where 
then matched to band 4. The radius of 1.5 pixels was 
chosen as a result of trying several different match radii 
between 0.5 and 2.5 pixels, and examining the number 
of matches between each band. If the match radius is 
smaller than the mean position uncertainty due to ran- 
dom errors, the number of matches increases rapidly with 
increasing match radius. This rate of increase falls off 
when the increase in the number of matches is mostly 
due to chance coincidences. The optimum radius is, of 
course, a function of band, but the similarity of the PSFs 
in the IRAC bands is enough that a single radius of 1.5 
pixels is a good pick for all four bands. The fluxes are 
"best" fluxes as described above with the aperture set to 
that of the shortest wavelength band in which the source 
is detected. The greater of the single-band catalog limit, 
or a three sigma limit were placed on non- detect ions (or 
negative aperture flux densities). The detection flag field 
in the four-band catalog has bit set for a channel 1 de- 
tection, bit 1 for channel 2, bit 2 for channel 3, and bit 
3 for channel 4. 

5. POSITIONAL ACCURACY 

For 2MASS sources in the main field detected in all four 
bands above the IRAC catalog flux limits (3223 sources; 
typical flux densities of » ImJy in channel 1), the mean 
radial position error with respect to 2MASS is 0.25 , 
with no measurable mean offset. This degrades at faint 
flux levels due to both noise and source confusion. Close 
to the flux limits of the catalogs our simulations suggest 
the positional error is w 1 . 

6. COMPLETENESS 

A series of simulations were run to estimate the com- 
pleteness of the main, verification, and ELAIS-Nl sur- 
veys. We define completeness as the chances of detect- 
ing a source placed at a random position in the field. 
Many sources are partially-resolved even at the limit of 
the surveys. To account for this in the completeness sim- 
ulations we needed to obtain an estimate of the range of 
source sizes near the survey limits. We used i-band Hub- 
ble Space Teleseope images from the Advanc;cd Camera 
for Surveys (ACS) taken as parallel data in program GO- 
9753 (P.I. Storrie-Lombardi) for which mosaics had been 
produced by I. Drozdovsky (p(TSonal communication) to 
study the morphologies of objects in the main and verifi- 
cation fields. For ELAIS-Nl no HST data were available. 



limits of the channel 1 and channel 4 catalogs of each 
of the three surveys (main, verification and ELAIS-Nl) 
were fit using software originally developed for quasar 
host galaxy fitting (Lacy et al. 2002). The results of this 
analysis are shown in Table 5. At the flux density lev- 
els of the main and verification surveys, stars contribute 
~ 30% of the sources, but these disappear in the ELAIS- 
XFLS softwarfj^]^ survey. The median galaxy half-light radius (ri/2) is 

XFLS software ," . , . , -n ■ r- n- r> i " 

jjgjjjj 0.4 m the mam and vcrificatKm surveys tailing to < 0.1 

in the channel 1 ELAIS-Nl data, and 0.25 in the chan- 
nel 4 ELAIS-Nl data. The distribution of galaxy types 
is field-galaxy like, with fa 20% ellipticals, « 80% spi- 
rals/irregulars. 

The scale sizes of sources in the main and verification 
fields mean that a significant fraction of these objects 
will be marginally resolved by Spitzer. In channel 1, 
the FWHM of the PSF is l.s", simflar to the FWHM 
of an exponential disk galaxy with ri/2 = 0.5 . In the 
ELAIS-Nl field, however, the scale sizes are sufficiently 
small that they can be treated as point sources for the 
purposes of completeness analysis. We therefore only 
modelled extended sources in the main and verification 
fields. For modelling the completeness in these fields, 
artificial sources were added to the mosaics with fluxes, 
sizes and proflles which are representative of the IRAC 
galaxy population close to the survey limit based on the 
averaged results of our fitting. Specifically, we picked 
the model population as follows: 30% point source, 60% 
i"i/2 = 0.4 disk galaxies and 10% ri/2 = 0.4 de Vau- 
couleurs galaxies. For channels 3 and 4 close to the main 
field flux limit, we increased the fraction of point sources 
to refiect the higher number of stars, using a mixture of 
60% point sources, 30% ri/2 = 0.4 disk galaxies and 

10% ri/2 = 0.4" de Vaucouleurs galaxies. We then es- 
timated the completeness by comparing the number of 
these artificial objects of a given flux density appearing 
in the Sextractor catalogs with the number known to 
have been added to the mosaics. The completeness was 
not a strong function of point versus extended source for 
the range of models used, or of model type. The com- 
pleteness plots for the main and veriflcation fields are 
shown in Figure 3. The ELAIS-Nl field has a very high 
source density in channels 1 and 2. For this field we 
adopted a slightly different strategy by picking a typi- 
cal field source from the mosaic, scaling it by different 
factors and inserting 1000 clones of this source into the 
mosaic. The results are shown in Figure 4. 

We also tried using the verification data to check the 
completeness of the main field catalog, by examining the 
fraction of veriflcation fleld sources which were detected 
in the overlapping region of the main fleld catalog as a 
function of flux density close to the survey limits in the 
main field. This gave completeness values higher by 10- 
15% close to our adopted survey limits. This result was 
not unexpected because a bright source will cause source 
confusion in both the main and the verification field, but 
provided a useful check on the simulations. 
We would like to emphasise that our completeness esti- 
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Fig. 3. — Results of completeness tests on the main and verification fields. The results for the main field are shown as open symbols, 
those for the verification field as closed symbols. The panels are for: (a) channel 1, (b) channel 2, (c) channel 3 and (d) channel 4. 



TABLE 5 

Morphological breakdown for XFLS sources close to the survey limits 



Channel flux range Fraction Median ri/2 range in ri/2 





(MJy) 


confused point disk de Vaucouleurs 


(arcsec) 


(arcsec) 


1 


20-30 


0/10 


2/10 6/10 


2/10 


0.4 


0.2-0.7 


1 


10-20 


0/10 


3/10 7/10 


0/10 


0.4 


0.3-1.1 


1 


3-4 


1/10 


0/10 8/10 


1/10 


0.07 


0.03-0.18 


4 


100-120 


0/10 


6/10 3/10 


1/10 


0.3 


0.2-1.4 


4 


30-40 


1/10 


3/10 5/10 


1/10 


0.4 


0.2-0.8 


4 


10-20 


1/10 


0/10 7/10 


2/10 


0.25 


0.1-0.5 



mates arc only approximate, particularly for the ELAIS- 
Nl channel 1 and 2 data, which is close to the confu- 
sion limit. A more careful analysis of completeness in 
this confused regime (see, e.g., Chary et al. 2004) will be 
necessary for accurate estimates of source counts at the 
faintest flux density levels in these data. 

7. PHOTOMETRIC ACCURACY 



An estimate of the uncertainty is associated with each 
flux density measurement in the catalog. This includes 
the uncertainties output by Sextractor, which are statis- 
tical in nature, and the known sources of systematic un- 
certainties. These have been quantified by the IRAC In- 
strument and Instrument Support teams (personal com- 
munications, though see discussions in the IRAC Data 
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Fig. 4. — Results of completeness tests on the ELAIS-Nl field. The panels are for: (a) channel 1, (b) channel 2, (c) channel 3 and (d) 
channel 4. 



Handbook), and are expected to amount to « 10%. Most 
of this is due to color-dependence in the flat field, and 
is manifest as scatter in the cataloged fluxes. This value 
of 10% was thus added in quadrature to the uncertain- 
ties from Sextractor. Wc have used three approaches to 
checking the systematic photometric errors. The recov- 
ery of fluxes from artificial sources was used to examine 
the uncertainties and biases in source measurement, in 
particular those associated with faint, slightly extended 
sources due to the small apertures and the consequent 
uncertainty in the aperture correction. Comparison was 
also made between fluxes measured for objects in the 
main field and the (much deeper) verification field. Fi- 
nally, comparison was made to 2MASS photometry of 
field stars. 

7.1. Fluxes of Artificial Sources 

The catalog fluxes of the artificial sources used in the 
completeness tests were compared to their true fluxes. 



This tests the accuracy and effectiveness of our aperture 
corrections and photometric algorithms, but does not ac- 
count for the color-dependent flat-field effect. 

These tests showed that, on average, more than 90% of 
the flux is recovered for point-like or small sources close 
to the flux limit of the catalogs. Below the flux limits, 
the fraction of the model flux recovered falls fast. The 
most likely reason for this bias is source confusion and 
variations in the background level of the mosaic affecting 
the local background estimate used by Sextractor. For 
recovered objects which were not blended with bright 
neighbors the difference between their true and recovered 
fluxes was in line with the expected statistical error. 

7.2. Comparison Between the Main and Verification 
Surveys 

Fluxes of objects in the part of the main survey which 
overlapped with the verification strip were compared. 
This comparison allows a rough estimate of the color- 
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dependent flat field error, as the main and verification 
surveys were taken with difi'erent grids, and of errors in- 
duced through over-zealous cosmic ray rejection in the 
shallower coverage main field data. 

The results of this test showed that the scatter between 
the main field and verification strip fiuxes is about 10% in 
all four bands for high signal-to-noise detections, roughly 
as expected from the known sources of systematic error. 
To get a true estimate of the error would require a more 
random superposition of the main and verification sur- 
vey grids, thus this estimate (which corresponds to an 
f» 7% uncertainty if the main and verification regions 
are affected similarly) is strictly a lower limit. 

The main field fiuxes in all bands average a little lower 
than those in the verification survey, by about 2%. This 
is probably due to a combination of factors. At the faint 
end, the fluxes are measured in small apertures in the 
shallower survey, frequently missing flux. At the bright 
end, the lower coverage of the main survey will lead to 
more frequent mis-identification of the brightest object 
pixels as cosmic rays, and the consequent exclusion of 
the brightest object pixels from the final mosaics. 

7.3. Colors of 2MASS Stars 

We performed a further check on the photometry by 
studying the colors of field 2MASS stars in the XFLS 
main field, following Eisenhardt et al. (2004). We placed 
a color cut oi J — K < 0.3 on the 2MASS stars, cor- 
responding approximately to the color of a F7 dwarf 
(Bessell & Brett 1988), and only considered stars with 
< 14 to ensure accurate 2MASS colors. Also, we 
excluded stars with < 10 from the quantitative com- 
parison in channel 1 as these were saturated in the IRAC 
images. For channels 2, 3 and 4 the mean stellar color is 
close to zero, as expected. In channel 1, our initial cal- 
ibration showed a significant offset from zero, similar to 
that noted by Eisenhardt et al.. Part of this is due to the 
expected non-zero colors of the stars. Most of the stars 
are close to the red limit in J — K, which corresponds 
to a K - L !v K - chl color of « 0.04 (see Table 6). 
However, a significant offset ~ 5%-8% still remained. 

Based on these plots, we decided to change the flux 
conversion factor for channel 1 to the value determined 
from measuring the fluxes of the A-star calibrators only 
(T. Megeath, personal communication), 0.1085, 4% dif- 
ferent from that given in the S9.5/S10.5 BCD headers 
(0.1125). This improves the channel 1 fluxes further, 
and the residual discrepancy is ~ 3%. This residual is 
comparable to that expected from the effects described 
above, and to the known uncertainties in the channel 1 
flat fields from IRAC campaigns 1 and 2. Figure 5 shows 
the final plots, and Table 6 the numerical results. Based 
on these, the fiux density scale of channel 1 looks accu- 
rate to « 3% both when compared to i^-band (panel (b)) 
and to channel 2 (panel (c)). The apparent upswing in 
the mean color in the K — ch2 plot for K > 13 (panel 
(d)) largely disappears when a redder color cut (which 
includes more objects) is used, so is probably not signifi- 
cant. Examination of Table 6 suggests that, indeed, our 
estimate of 10% systematic uncertainty in the flux densi- 
ties of individual sources is reasonable, and that the sys- 
tematic uncertainty in the flux density scale is ~ 2 — 3%. 
Unfortunately there were too few stars in the verification 
region to allow us to perform the same test on the deeper 



TABLE 6 

Offsets and dispersion of the colors of blue stars 

BETWEEN 2MASS AND IRAC 

Color Expected Offset from Dispersion around 

mean value* expected mean expected mean 

K - chl 004 +003 O07 

K-ch2 0.02 -1-0.00 0.12 

K-ch3 0.00 -0.02 0.11 

K-chi 0.00 -1-0.02 0.12 
• from Bessell and Brett (1988). 



data, but the comparison between the main and verifi- 
cation field fiux densities showed only small systematic 
differences (Section 7.2). 

8. SETTING THE FLUX DENSITY LIMITS 

As this survey is expected to be used primarily for sta- 
tistical purposes, we have produced single-band catalogs 
cut at flux density levels at which the surveys are still 
« 80% complete (Table 7), and at which the simulations 
indicate that the mean fraction of flux recovered by Sex- 
tractor from point or slightly extended sources is 90%. 
This resulted in us setting limits in the main field cata- 
log of 20/iJy, 25/iJy, 100/iJy and 100/iJy in channels 1, 

2, 3 and 4, respectively. The corresponding limits in the 
verification survey were 10/iJy, 10/iJy, 30/LtJy and 30/zJy, 
and those in the ELAIS-Nl field are 4/iJy, 3/xJy, 10/zJy 
and 10/iJy. The survey limits scale approximately with 
the square root of the exposure time in channels 2-4, 
but channel 1 is less deep than expected in the verifica- 
tion and ELAIS-Nl fields. Source confusion is certainly 
partly to blame, but in channel 2 similar source den- 
sities are reached in the deep surveys without a strong 
departure from the expected scaling. The most likely ex- 
planation is a combination of source confusion and the 
latent image problem, which was especially severe in the 
verification data in channel 1, and resulted in a vary- 
ing background which was only partly corrected by the 
subtraction of a delta-dark. The fiux density limits cor- 
respond to a signal-to-noise of « 5 for a typical source at 
the limit of the surveys (Figure 6), except in the deeper 
channel 1 data affected by variable background and con- 
fusion, where it is « 7. 

We have compared our fiux density limits with those 
of the Spitzer Performance Estimation Tool (PET, 
http://ssc.spitzer.caltech.edu/tools/pet.html). For the 
main field, we expect ba sensitivities of 7.4, 11, 61 and 
60/zJy in channels 1-4 respectively, compared to our sur- 
vey limits of 20, 25, 100 and 100/iJy, respectively. In the 
verification survey, the PET gives 5a sensitivities of 1.8, 

3, 20 and 23 /iJy compared to 10,10,30 and 30 /iJy, and 
in ELAIS-Nl, 0.33, 0.7, 4.6 and 5.8 ^J.3y compared with 4, 
3,10 and 10/iJy. There are three reasons why we do not 
achieve the PET sensitivities. First, the PET estimates 
are based on point-source fitting rather than aperture 
photometry. Point-source fitting has 7-13 noise pixels, 
depending on channel (see the Spitzer Observers' Man- 
ual), compared to 19 pixels in our smallest photometric 
aperture, which accounts for much of the difference in the 
main field survey, and in all surveys in channels 3 and 

4, Second, the PET estimates take no account of source 
confusion. This is particularly important for channels 1 
and 2 in the ELAIS-Nl survey. Third, the background 



10 



Lacy et al. 



variations in channel 1 discussed above limit our sensi- 
tivity in the deeper surveys in that channel. 

9. RELIABILITY 

Objects detected in the main field catalogs which over- 
lapped the region covered by the verification data were 
compared with those made from the 30s data in the ver- 
ification strip, which goes much deeper. This allowed us 
to make an estimate of the reliability of the main field 
catalog at the catalog flux density limits given above. 
These tests showed that the catalogs are « 99% reliable 
in all four bands, in the sense that the probability of an 
object in the main field catalog being detected in the 
deeper verification data was ~ 99%. Measuring the reli- 
ability of the verification and ELAIS-Nl catalogs is not 
possible using this technique, but as their limits are sim- 
ilar to the main field in terms of signal-to-noise ratio we 
expect these catalogs to be similarly reliable. 

10. DATA PRODUCTS 

10.1. Mosaics 

A mosaic for each channel for the main and verifica- 
tion fields is available, along with coverage maps. All 
are made in the same fiducial frame, i.e. with the same 
pixel grid. This results in much of the verification mosaic 
being blank, but which greatly eases source comparison 
as the same pixels in every mosaic and in every channel 
correspond to the same sky position. The 12s and 30s 
frames in the verification region are coadded according to 
inverse variance weights determined from the measured 
standard deviations in the BCDs. We opted to keep the 
main field data separate from the verification field data 
to maintain their statistical independence. The coverage 
maps show the numbers of frames added to produce the 
final mosaic. In the main field, one unit of coverage cor- 
responds to one 12s frame, in the verification field one 
unit is one 30s frame, with the 12s frames being 0.388 
(the ratio of the exposure times in the 12s and 30s data). 
Multiplying by the appropriate exposure time (10.4s for 
the 12s frames in the mainfield and 26.4s for the 30s 
frames in the verification field) will thus recover the ex- 
posure time for a given pixel. 

10.2. Catalogs 

10.2.1. Single band catalogs 

The single-band catalogs contain the following fields: 
srcID - an integer identifier for the source, x, y - 
pixel coordinates of the source in the mosaic, RA, 
Dec - right ascension and declination of the source in 
the single-channel mosaic (J2000 coordinates), best flux 
most reliable estimate of the flux density selected as 
described in Section 4, best aper - the aperture for 
the "best" flux (in arcsec), err best flux - the uncer- 
tainty in the "best" flux, S:N - the signal-to-noise ra- 
tio in the "best" aperture, 6'.'00 ape - aperture-corrected 
flux density in a 6'.'00 diameter aperture, 9'.'26 ape - 
aperture-corrected flux density in a 9'.' 26 diameter aper- 
ture, 14'.'86 ape - aperture-corrected flux density in a 
14'.'86 diameter aperture, 24!.' 4 ape - flux density in a 
24'.'4 diameter aperture (no correction applied as this is 
the aperture the IRAC calibration stars are measured 
in), err G'.'OO, err 9'.'26, err 14'.'86, err 24'.'4 - uncertainties 
in the aperture flux densities, area - area of the isophotal 
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Completeness and confusion statistics 
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TABLE 8 

Numbers of sources in each catalogue as a function of 



detection flag 
binary (decimal) 


main 


BAND 

verification ELAIS-Nl 


0001 (1) 


30758 


2631 


1271 


0010 (2) 


6984 


2061 


1622 


0011 (3) 


49534 


3610 


106 


0100 (4) 


822 


139 


655 


0101 (5) 


698 


666 


627 


0110 (6) 


24 


16 


2 


0111 (7) 


2464 


582 


34 


1000 (8) 


821 


114 


297 


1001 (9) 


61 


8 


1 


1010 (10) 


727 


525 


543 


1011 (11) 


3510 


368 


25 


1100 (12) 


103 


12 


3 


1101 (13) 


11 


2 


4 


1110 (14) 


23 


18 


4 


1111 (15) 


6653 


1472 


45 


total in catalog 


103193 


12224 


5239 



aperture (in pixels), Ixx, lyy, Ixy- second moments of the 
flux distribution, flag - sum of the flag bits, see Table 3, 
class - stellarity index from Sextraetor (between and 1, 
low values correspond to probable extended sources and 
high values to point-like sources), coverage - number of 
unit frames contributing to the mosaic at the source po- 
sition (see section 10.1). All fluxes and errors are quoted 
in /iJy. The errors include the statistical error and the 
10% systematic uncertainty discussed in Section 7. The 
uncertainty in the overall flux density scale is not in- 
cluded, but we believe this to be small in any case, see 
section 10.3). Table 7 summarizes the properties of the 
single-band catalogs. 

10.2.2. Four-band catalogs 

The four-band catalogs contain the following fields: 

x,y - pixel coordinates of the source in the shortest- 
wavelength mosaic in which it was detected, RA, Dec 
- right ascension and declination of the source in the 
shortest-wavelength mosaic in which it was detected 
(J2000 coordinates), chl id, ch2 id, ch3 id, ch4 id - iden- 
tifiers in the four single-band catalogs (defaults to -1 for 
non-detections in a given band), chl flux, err chl - chan- 
nel 1 flux density and uncertainty in the "best" aperture. 
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Fig. 5. — Photometry of blue stars in the XFLS field, (a) J—K versus K color-magnitude diagram. The dashed line shows the J—K = 0.3 
color cut used to select the blue stars (corresponding to late-F and bluer), (b) K — chl vs K. The dashed line at K — chl = 0.04 is the 

expected color of stars with J — K = 0.3. (Stars brighter than K ^ 10 are saturated in the channel 1 image, producing the downturn at 
bright magnitudes.) (c) chl — ch2 vs K. The dashed line at chl — ch2 = —0.02 is the expected color of stars with J — K = 0.3. (d) K — ch2 
vs K. The dashed line at K — ch2 = 0.02 is the expected color of stars with J — K = 0.3 , {e) K — ch3 vs K, and (f) K — chA vs K. In (e) 
and (f) the dashed line is at zero, the expected color in these pairs of bands. 



ch2 flux, ch2 err channel 2 flux density and uncertainty, 
chS flux, chS err channel 3 flux density and uncertainty, 
ch4 flux, ch4 err channel 4 flux density and uncertainty, 
aper - aperture used for flux density measurements, flag 
- see Table 4, dflg - sum of the detection flag bits, see 
discussion in Section 4. 

Objects with no coverage in a given band have their 
flux densities in that band set to 99999.0. Otherwise, 
non-detections have their fluxes set to the appropriate 
catalog limit, or a '3-a limit, whichever is higher. 

Table 8 shows the numbers of sources in each catalog 
as a function of detection flag. Of 93689 channel 1 main 
field detections, 62161 (66%) are also detected in channel 
2, 9826 (10%) arc detected in channel 3 and 10235 (11%) 
are detected in channel 4. 

11. LESSONS LEARNED AND FUTURE WORK 

The observations for the XFLS were designed well be- 
fore launch, when the knowledge of the properties of the 
instrument was rather limited. With the benefit of hind- 
sight we would have changed a number of aspects of the 
design. The greater than expected scattered light prob- 
lems, and the latent image problems in channel 1 both 
could have been mitigated by using a larger dither pat- 
tern (with arcminutc-scalc dithers) and/or half array off- 
sets in the mapping strategy. On the other hand, one 
correct decision was to maintain a relatively high cover- 
age factor (five over the main field). This allowed reliable 
cosmic ray rejection, and the extra redundancy allowed 
us to use the scattered light masks effectively. 

Future reprocessings of the XFLS data should be able 



to improve the accuracy of the fluxes. For example, we 
may be able to use color-dependent flat fields to obtain 
more accurate fiuxes by applying the flat field appropri- 
ate to the color of each object using an iterative tech- 
nique. Also, we may be able to use improved source ex- 
traction techniques based on a development of the SSC 
point source extractor to optimally-extract the typically 
slightly extended IRAC sources, and properly quantify 
the confusion in the ELAIS-Nl channel 1 and 2 data. 
These improvements were considered to be beyond the 
scope of this work, where our purpose is to provide a re- 
liable and fairly complete catalog on a reasonably short 
timescale. 

12. SUMMARY 

The IRAC data from the XFLS has been analyzed, and 
a set of catalogs produced which we beheve to be « 80% 
complete and « 99% reliable. The final bandmerged cat- 
alogs contain 103193 objects in the main field, 12224 in 
the verification field and 5239 in ELAIS-Nl. Flux densi- 
ties of high signal-to-noise objects are accurate to about 
10%, and the systematic uncertainty in the absolute flux 
density scale is 2 — 3%. Positional accuracy is w 0'.'25 
for high signal-to-noise objects and « 1 at the flux den- 
sity limits of the catalogs. We have successfully extracted 
sources at source densities as high as 100000 deg^^ in 
our deepest channel 1 and 2 data, though there are in- 
dications that we are approaching the confusion limit at 
these high source densities, in agreement with Fazio et al. 
(2004b). The mosaics and catalogs will be made avail- 
able both through the Spitzer Science Archive and the 
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Fig. 6. — Flux density against signal-to-noise ratio in the "best" apertures for sources in the raw main field catalogs in each channel: (a) 
channel 1, (b) channel 2, (c) channel 3, and (d) channel 4. The dashed lines show the adopted fiux density limits for the final catalogs in 
each band. The striations in the plots are the result of selecting the fluxes from a small set of fixed apertures. 
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